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Abstract. 

It has been recently shown that collisional models can be used to derive a general 
form for the master equations which describe the reduced time evolution of a composite 
multipartite quantum system, whose components "propagate" in an environmental 
medium which induces correlations among them via a cascade mechanism. Here we 
analyze the fundamental assumptions of this approach showing how some of them can 
be lifted when passing into a proper interaction picture representation. 
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1. Introduction 

In the study of the open dynamics of a multipartite quantum system it is frequently 
made the simplifying assumption that each subsystem interacts with its own local 
environment. In the language of quantum communication [T| this is equivalent to 
saying that the resulting time evolution is memoryless i.e. that the noise tampering 
the communication acts independently on each local component (information carrier) 
of the transmitted quantum message. In recent years, however, it has been shown that 
interesting new features emerge when one makes the realistic assumption that the action 
of a channel over consecutive uses is correlated [21 El IH El El El IS]- Such correlations 
have been phenomenologically described in terms of a Markov chain which gives the joint 
probability distribution of the local Kraus operators acting on the individual carriers [2] . 
Alternatively they have been effectively represented in terms of local interactions of 
the carriers with a common multipartite environment which is originally prepared into 
a correlated (possibly entangled) initial state [H], or with a structured environment 
composed by local and global components O HI E] • These models, although physically 
well motivated do not have an intrinsic time structure, in other words they are unable 
to describe a situation in which the information carriers interact one after the other 
with an environment which in the meanwhile evolves. For instance consider the case in 
which an ordered sequence of spatially separated photon pulses carrying information in 
their photon number, propagates at constant speed in a non passive, lossy optical fiber 
characterized by (relatively) slow reaction times. If the speed of the pulses is sufficiently 
high, one might expect that, thanks to the mediation of the fiber, excitations from one 
pulse could be passed to the next one modifying their internal states via a (partially 
incoherent) cascade mechanism (see Fig. [l] a) for a schematic representation of the 
process). The net result of course is the creation of delocalized excitations over the 
whole string of carriers as they proceed along the fiber (the extent of such delocalization 
depending upon the ratio between the transmission rate at which the carriers are fed 
into the fiber and the dissipation rate of the latter). Alternatively, consider the case 
in which a linear array of local quantum systems (say a set of driven QED cavities as 
in Fig. [l] b), or a set of quantum dots composing a quantum cascade laser [U]), are 
indirectly coupled via unidirectional environmental mediators (the photons emitted by 
the cavities or by the dots) which passing from one system to the other, allow them to 
exchange excitations [TUl HH [121 [131 [IH [ISl [IS]- As in the previous case, the formation 
of delocalized excitations is expected as time passes (in this case the delocalization of 
the excitation will depend upon the product between the damping rate of the mediator 
and the distance between two consecutive quantum systems). 

The general form for master equations which describe these situations has been 
recently derived in Ref. [17] by adopting a collisional approach [TSl HH] to describe 
the system/environment coupling. We aim to review these findings, focusing on some 
technical aspects of the problem which allows us to lift some of the assumptions of 
Ref. p7]. The manuscript is organized as follows. In Sec. |2]we describe the collisional 
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Figure 1. (Color online) Prototypical examples of quantum systems which admit 
a representation in terms of cascade quantum master equations, a) A string of 
quantum information carriers (say the photon pulses represented by the green element 
of the figure) propagating along a fiber while losing photons along the way. Here 
the sub-environment Ei,E2, - ■ ■ are associated with different sectors of the fiber which 
effectively describe the transmission line within a lumped model approach. If the 
propagation speed of the pulses is sufficiently high, the excitations they waste in a given 
sector of the fiber cannot be re-absorbed by the same pulse: still the subsequent pulses 
have a certain probability of absorbing it. The net result is an indirect, unidirectional 
coupling between the pulses mediated by the fiber which allow quantum signals to pass 
from one carrier to the subsequent ones as schematically shown by the yellow arrows 
of the figure, b) A set of atoms or ions (represented by the orange elements in the 
figure) trapped into a series of QED cavities which exchange photons (blue arrows), 
via unidirectional couplers (black elements) [lOl [11] . c) An array of cavities , S'2 , • • • 
crossed sequentially by atoms Ei,E2, - ■ ■ (blue elements) of a beam. The injection rate 
is such that the atoms cross the cavities one by one. The atoms are initially prepared 
all in the same state. 



model, its continuous time limit (Sec. 2.1) and the basics properties of the associated 
master equation for cascade quantum systems. In Sec. [3] we then pass to discuss 
the fundamental assumption which underline the derivation (namely the environment 
stability condition under collisions). Here we first show how local free evolution term 
can be embedded in the derivation (Sec. 3.1 ). Then we prove that the stability condition 
can always be enforced by passing through an interaction picture representation which 
defines a more "stable" effective coupling with the environment (Sec. 3.2). Conclusions 
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Figure 2. (Color online) Schematic of the collisional model. The horizontal lines 
describe an ordered set of carriers 5i , 5*2 , • • • which interact with an ordered set 
of (possibly infinite) identical local sub-environments Ei, E2, ■ ■ ■ via local unitaries 
Us^E„ represented by the yellow elements [r] being the initial state of the Ejs). 
Between collisions each sub-environment evolves according to the "damper" maps M. 
(represented in the figure by the grey elements). The overall dynamics can be described 
as a ordered sequence of row or of column super-operators - see Fig. [3] for details. 



2. The collisional model 

Consider a multipartite quantum system 5, composed by M - not necessarily identical 
- ordered subsystem Si, S2 ■ ■ ■ Sm which we shall refer to as the information carriers 
of the model. They are assumed to be prepared in a possibly entangled initial state 
p(0) and to evolve in time due to the interactions with a multipartite environment 
S consisting of a collection of sub-environments Ei, E2, ■ ■ ■. Following the collision 
model of irreversible dynamics presented in Refs. flSl fT9]. the carriers/sub-enviroment 
couplings are described via a sequence of pairwise, time-ordered unitary interactions 
which in our case are organized as in the scheme shown in Fig. |2} According to it, each 
element of S interacts with all the elements of £ in such a way that given m' ^ m 
and n' ^ n integers, the "collision" between Sm and En is assumed to happen before 
the one involving Sm' and En'- In particular, this implies that the Si,Ei interaction 
takes place before the couplings between Si,E2, the coupling between 5*2, and the 
coupling between 5*2, i?2- Similarly the 5*2, £'2 coupling is assumed to come after the 
5*2, El and the Si, E2 couplings, while no specific ordering is imposed on these last two 
events. Within this theoretical framework the temporal evolution of m-th carrier Sm 



and final remarks are presented in Sec. |4j 
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can be then be described trough the action of the following joint unitary evolution 

where, for instance, 

Us^,E,. := exp[-tgHs^E^At] , (2) 

is the transformation that characterizes the "collision" between 5*^ and Em- In this 
expression At is the collision time, g is an intensity parameter that gauges the strength of 
the coupling, while Hs^En is the coupling Hamiltonian which, without loss of generahty, 
we write as 

J^S„.,.:=E'»"«4',:. (3) 

e 

with Ag\B]^^ being Hermitian operators. It is worth stressing that in writing Eq. 
one implicitly assumes that a given carrier never interact twice with the same sub- 
enviroment. This hypothesis is typically enforced in collisional models which aim to 
describe Markovian processes - see however Ref. [20j for an alternative approach. Its 
validity relies on the existence of a two well separated time scales: a fast one, which 
defines the typical correlation times of the environment, and a slow one, which instead 
defines the dissipative effects on the system of interest (i.e. the carriers) induced by 
the coupling with the bath. Such assumption of course is not always fulfilled and 
when enforced it inhibits the possibility of feedback mechanisms where the state of the 
system at a given time is influenced by the entire evolution history. Notice however 
that in the scenario we are considering here, the global Markovian structure of the 
coupling ([T]) doesn't prevent the possibility that different carriers could have a non 
trivial causal influence on each other. In other words, as schematically shown in Fig. [T] 
a), quantum "information" can be transferred from one carrier to the other through the 
intermediation of the environment. 

In top of the processes described by the unitary couplings ([T]) we also assume that 
between two consecutive collisions each sub-environment evolves according to the action 
of a local Completely Positive (CP) map Ai. The latter is introduced to effectively 
account for the internal dynamics of in particular the transformations Ai mimics the 
relaxation processes that may take place within the environment alone (e.g. originating 
from the mutual interactions between its various parts) and which in principle involve 
timescales different from those that define the rate of the collisional events. In other 
words, as in Ref. |5], the mappings Ai act as "dampers" for the information that 
percolates from one carrier to the subsequent on^ how effective such damping is, it 
depends of course upon the rate at which two subsequent carriers approach the same sub- 
environemnt (i.e. in the example of Fig. [l]a), it depends upon the propagation velocity 
of the pulses along the fiber). Putting all together the resulting temporal evolution can 

I In what follows we will work under the simplifying assumption that the same CP transformation 
acts among any two collisions - the generalization to the case in which the A4 change passing from one 
collisional event to the other being straightforward. 
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Figure 3. (Color online) Schematic of the decomposition of the evolution map 
collisional model in terms of a) row operators as detailed in Eq. ([5| (in the figure these 
operators are represented by the gray elements), or b) column operators (represented 
by the red elements) as detailed in Eq. Q. 



hence be expressed in a compact form by observing that after the interactions with the 
first n elements of £ the global state R{n) of the system and of the environment is 
obtained from the initial state p(0) ® rj®^ as 

/?(n) = >V("'*')(p(0)®r/®") , (4) 

where W^"'*^) is the super-operator which describes the collisions and the free evolutions 
of 8 while rj is the density matrix which describes the initial state of the sub- 
environments (for simplicity we assumed that all the En are characterized by the same 
initial state). The VV*^"'^-' can be expressed as a composition of row super-operators 
stack in series one on top of the other (see Fig. |3]) 

y^iuM) = , o • • • o Teg, o Teg, , (5) 

where we use the symbol "o" to represent the composition of super-operators and where 

Tef) o Wg, . (6) 

In the above expression, given a unitary transformation U , we define Ll{- ■ ■) = U{- ■ ■)U\ 
while we used the symbol A^*^" to represent A^^;^ o ■ ■ ■ o A^^^, -^Ej being the map Ai 
operating on the j-th element Ej of £. The transformation , describes the evolution 
of Sm in its interaction with S plus the subsequent free evolution of the latter induced 
by the maps Ai. Alternatively, exploiting the fact that for m' m, n' n the operators 
Usm,E„ and Us^,,e^, commute, W'^"'*^^ can also be expressed in terms of column super- 
operators concatenated in series as follows: 

>V("'"^=C^oC^_^o... 045)04^, (7) 
where for all j = 1, ■ ■ ■ , ra, 

■■= Me, o Us,„e, o...oMe,o Us„e, ■ (8) 
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Equations ([s]) and ([s]) enlighten the causal structure of the model. In particular ^ 
makes it explicit that whatever happens to 5*^+1 comes always after the transformations 
operating on Sm- As a consequence the latter can have indirect influence on the former 
but the vice-versa is not allowed. Similarly Eq. ([s]) shows that an analogous causal 
structure is present on the elements of £: events involving En may have causal influence 
on those involving En+i but the opposite is impossible. This last equation is also useful 
to write a recursive expression for R{n). Indeed by construction we have 

i?(n+l)=4^)^^^(/?H®r^), (9) 

which confirms the intrinsic Markovian structure of the temporal evolution for the whole 
set of carriers that composes S. The recursive form of Eq. (|9]) is similar to the one 
characterizing the models of Refs. [19] where, for a single qubit carrier (M = 1) and for 
a particular class of interaction unitaries, it was shown that it leads to a dynamics which 
can be described by a Lindblad super-operator. Following Ref. [IT] one can generalize 
this fact to an arbitrary number of carriers and for arbitrary coupling Hamiltonians ([S]). 
We simply assume a weak coupling regime where we take a proper expansion with 
respect to the parameters g and At which quantifies the intensity and the duration of 
the single events. In particular we work in the regime in which g^t is small enough to 
allow for the expansion of the dynamical equation ^ up to (9(((y'At)^), i.e. 

R{n + l) = [Xs,E^^,+C's,E^^^g^t (10) 
+ C'^^E^^^{g/\tf]{R{n) ®r^) + 0{[g^tf) , 

where Ts,e„+i is the identity superoperator while C^^^^^^ and C^^^^^^ are the first and 
second expansion terms in gAt of the superoperator C^^^^^_^, respectively (see below). 
The resulting expression can then be traced over the degree of freedom of S to get an 
equivalent expression for the the reduced density matrix of S alone, yielding. 

Pin + 1) = Pin) + igAt) (i?(n) ® r^)>^ (11) 

+ igAtf {Clj,^jRin)®v)), + OUgAtf) , 

where we used the symbol (■ ■ ■)£ to represent the partial trace over Ei,E2, ■ ■ ■ and where 
for all n we introduced 

Pin) := {Rin))e , (12) 

(it represents the joint state of the carriers after the interaction with the first n sub- 
environment). Explicit expressions can be obtained by noticing that for each m and j, 
the super-operators Us^^Ej admit the following expansion, 

Us^,E, = Xs^,E, + igAt) U' + igAt)' W' + OUgAt)') , (13) 



with 



—I 



HSm,E 



m 5 5 



Wc' ■ ■) := Hq_,eA- ■ ■)Hq_,e, 



Jj2 

2 Sm^Ej ' 



(14) 
(15) 
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where [■■■,■■■]_ and [■■■,■■■]+ represent the commutator and the anti-commutator 
brackets respectively. From Eq. Q it then follows that 

M 

^s,E, ■■= E ^E^^' ° ^Le, ° M'^-' , (16) 



with 



M 

<m— 1 



'E, ' 



Cs,E, ■= 2^ ^E, ° ^S^,E, ° 

m=l 

M M-l 

^s%, ■■=1111 M%-'^' o U's^,^^^ o M%- o W^^^^^. o , (18) 

m'=m+l m=l 

(here stands for the iterated application of m maps on the same sub- 

environmental system E, e.g. A^|; = J^e ° -Me)- Replacing these expressions into 
Eq. (11) and remembering the definition (|3| the first order term in gAt gives 

{Cs,E„,M^)^v)), = -^[llH(:ff\p{n)]_ , (19) 

m 

with Hm^^'^ being the following effective local Hamiltonians 

HI-^^'--=11{bZ.Me:IM)^ 41- (20) 



For the second order terms in g/S.t we get instead two contributions associated 
respectively to local Lindblad terms (i.e. Lindblad terms which act locally on the m-th 
carrier) and two-body non local terms which couple the m carrier to the m' > m. More 
precisely the first one is given by 



^(i?(n)®r^)) =- J2^Mn)), (21) 



where 7 is a positive parameter whose value will be specified later (see Eq. (29) below), 
while Cm is the super-operator 



cu---) = llliii''W^li---)A^ 



- A^p{... )-{... )Af A^P]. (22) 
In this expression the coefficients 

1^1'''--=i{b'e'b^PmT\v))e, (23) 

define the (non negative) correlation matrix of the sub-environment operators and 
evaluated on the density matrix Ai''^^^{r]) which describes the state of the sub- 
environment after m — 1 free (i.e. non collisional) evolution steps. Equation (22) can 
also be casted in a more traditional form |2T] by diagonalizing 7^' • this allows one to 
identify the decay rates of the system with the non-negative eigenvalues Vm of 7m'^ "* and 
the associated Lindblad operators L^^ with a proper linear combinations of the A^f^ . 
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by 



The second contribution of order two in (7 At which enters Eq. (11 ) is instead given 



e 7 



(24) 



in'>m 



where for m' > m V^I^K is the super-operator defined as 



^m,m'\ J ~ 2-^ fm,m' ^S„ 
/ J L ' m,m'i 



in 



A 



{(■) 



(25) 



with 



7^2 ■■= 7 {B^pM-'--{B^^ M-^-\v)) )e . (26) 

The coefficients 7.)^'2' introduce cross correlations among the carriers and depend upon 
the distance m' — m between the associated rows of the graph of Fig. [2j Furthermore, 
similarly to the the terms of Eq. ( 23 ) , they also depend on m — 1 due to the fact that 
the model admits a first carrier. 

The resulting expression for p{n + 1) can thus be written as 

p(n + 1) - p{n) 



+ 



At 
^?'At 



7 



m 

J2^m{p{n)) + J2 '^i^liPin)) \ + 0{g'At') , (27) 



m'>m 



which possesses an explicit Markovian structure characterized by the presence of an 
effective Hamiltonian (first line) and dissipative contributions (second line). 



Before proceeding to the continuos limit let us briefly review how the above scheme 
applies to a specific discrete system, with the aim to clarify the meaning and the validity 
of the assumptions made in our model. For this purpose we refer to the prototypical 
example of Fig. [l] c). Here an array of cavities is driven by a beam of atoms crossing 
them. The rate of injection of the atoms is such that the atoms enter the cavity one by 
one as shown in figure. The atoms in the beam, all initially prepared in the same state 
cross sequentially all the cavities of the array. The atom-field interaction is described 
by the Jaynes-Cummings hamiltonian, which can be straightforwardly cast in the form 
(§. For a single cavity crossed by a beam of single atoms, in the absence of leackage of 
photons out of the cavity, such model describes damping or amplification for the cavity 
field, whose dynamics can be described by a Markovian master equation in Lindblad 
form [22] . The extension to n cavity is described naturally in our model. Such markovian 
behavior is due to two crucial features of the model: the short (finite) time At which 
takes each atom to cross a cavity and the fact that the atomic state is "refreshed" when 
a new atom is injected. The master equation so obtained describes a coarse grained 
time derivative on a timescale At. On such time scale the environment is reset. In the 
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standard theory of damping of a system which is continuosly interacting with the same 
- big - reservoir the time At would be the self correlation time of the reservoir. This 
is not the case in our scenario: each sub environment is small but it interacts for a 
short time and then, after a time At substituted with a new one. Furthermore the cross 
terms in our master equation do not describe a collective, simultaneous, coupling of the 
subsystems with the environment, but rather they describe how the dynamics of the 
various subsystems are correlated due to the fact that they have interacted sequentially 
with the sub environments. This explains why the dynamics described by our system 
is markovian and does not exhibit the non markovian multipartite features which are 
characteristic of the scenarios analyzed in [23] . 



2.1. Continuous limit 



Equation (27) can be turned into a continuos time expression by taking the proper limit 



At — > while sending n to infinity so that 

lim n At = t . (28) 

Notice that there are two possible regimes. If g is kept constant as At goes to zero. 



then the dissipative contributions of Eq. (27) are washed away and the dynamics 



reduces to a unitary evolution characterized by the effective (possibly time-dependent) 



Hamiltonian (20). The situation becomes more interesting if instead g is sent to infinity 



so that g'^At remains finite, i.e. [IT] 

2 



7 • (29) 
Enforcing this limit is of course problematic due to the presence of the first order 



lim g^At 

Ai-5-0+ 



contribution in Eq. (27) which tends to explode. A way out is to assume the following 



stability condition for the environment [17j . 



(5^ A^"^(r/))s = V£,m 



(30) 



which ensures that Hm'^^^ 



and hence the first order contribution of Eq. (27), identically 
nullifies. Under this hypothesis, defining p(t) = liim.At^o+ piji)i one can indeed arrive 
to the following continuous master equation for the system, 

p(t) = cMt)) + J2 ^ti'(pW) ' (31) 



m=l 



m'>m 



whose properties have been characterized in Ref. [T7]. Here we only mention that the 



cross terms appearing in Eq. (31 ) have an intrinsic unidirectional character which makes 



this expression suitable to characterize the dynamics of a cascade quantum system. 



Indeed, for each m' > m it can be directly verified from Eq. (25) that we have 

vt'li---)) =0. 



(32) 



This implies that the evolution of the first m carriers of the system are not influenced 
by the evolution of the ones that follow (in other words it is possible to write a master 
equation for the density matrix of the first m elements of S only). The opposite 
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relation however is not true as in general "P^^.^/ doesn't nullify when traced over Sm, 
i.e. {T^^mrn'i' ' ')/ 7^ ^- This means in particular that in our model it is in general 
impossible to write a master equation that involves only the density matrix of the m'- 
th carrier 5*^' (we need indeed to include also the carriers that precede it[|]). Similar 
properties were obtained in the works of Gardiner, Parkins and Carmichael (TU] in theirs 
seminal study of cascade optical quantum systems. As shown in Ref. [T7j the latter can 
be seen as special instances of (31) for specific choices of the couplings ^ and of the 
environment initial state rj. 



In what follows we will not discuss further the implications of Eq. (32). Instead 



we will focus on the assumption Eq. (30) showing how it can be enforced by passing in 



a proper interaction picture with respect to the free evolution of the carriers. Before 
doing so however we think it is worth stressing that the above derivation still holds also 
if the collisional Hamiltonians ([3]) are not uniform. For instance suppose we have 

(33) 



SmEn 



B 



En 



where now the operators acting on the carrier Sm are allowed to explicitly depends 
upon the n index which label the collisional events, and similarly the operators 
acting on the sub-enviroment are allowed to explicitly depends upon the index m 
which labels the carriers. Under these conditions one can verify that Eq. (27) 



is still valid even though both the super-operators and T>^1^, become explicit 
functions of the carriers labels and of the index n which plays the role of a temporal 
parameter for the reduced density matrix p{n). Specifically they are now defined 

instead of and 



respectively as in Eqs. (21) and ([24) with the operators A^^ 
with the coefficients {B^^P M"'-\r]))E and (sj^'^A^™'-™^ ^^^^ 
by (^Sr-'^^iT'''^ M^-\ii))e and 



(n+l,^) 



IE anu (5^'"''^'^7W™'-™(S{f'^) M'^-^{r,)))E respectively. 
Similarly the continuous limit can be enforced as before: in this case however, to account 



{B^^ replaced 



for the non uniformity of the couplings, the condition (|30|) becomes 



B 



0, Vm,£. 



Furthermore, while taking the limit (28) the operators A 



(34) 

acquire an explicit 

temporal dependence which transforms them into a one parameter family of operators. 
As a result we get a time-dependent master equation characterized by a Lindblad 
generators which explicitly depends on t, i.e. 

p{t) = j2^m{pm) + Yl T^mX'ipm) , (35) 



m=l 



m'>m 



with Cr 



t) and V'-^^ 



m,m' ^ 



replaced by Af (t) := limAt^o+ ^ 



(n+i,e 



■,t) as in Eqs. (22) and (25) with the operators A 



§ A notable exceptions being the case in which all the coefficients 7^'^'! appearing in Eq. (25) are real: 
under this condition (t>1^1^, (• • = so that the dynamics of every carrier is causally disconnected 

from the others [17] . 
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3. The stability condition 



The name stability condition given to the constraint (30) follows from the fact that 
it implicitly assumes that during the collisions the sub-environments are not affected 
by the coupling with the carriers (at least at first order in the coupling strength). 
This is mathematically equivalent to the standard derivation of a Markovian master 
equation [21] for a system interacting with a large environment, in which one assumes 
that the overall system-environment density operator at any given time t of the evolution 
factorizes as in p{t) ® rj where rj is the environment density operator. The two scenarios 
are however different. In the standard case the reason for which the environment state 
is unchanged is because it is big. In the scenario analyzed here, consistently with the 
collisional model, the environment state is constant because each subsystem collides 
briefly with a sequence of sub-environments all initially in the same state. 



As anticipated in the previous section, in our analysis of the condition (30) a 
proper handling of the carriers free evolution plays a fundamental role. This should 
not come as a surprise: an important step in the standard derivation of a Lindblad form 
is the possibility of effectively "removing" the free evolution of the system and of the 
environment by passing in the associated interaction representations. Such step is useful 
because it allows one to directly relate the fast evolution times of the large environment 
with the slow decaying rates of the system of interests: it is in this limit that the Markov 
approximation can be properly enforced In our model we can show that the cases 



in which Eq. (30) cannot be directly enforced, can be mapped into effective models in 
which Eq. (30) exactly holds but which allows for explicit free evolution terms for the 
carriers between any two collisions which have to be removed by passing in a proper 
interaction picture representation. As a preliminary step toward the discussion of the 
stability condition it is hence important to discuss how the derivation changes in this 
last circumstance. 



3.1. Including local free evolution terms for the carriers 



Assume that the stability condition (30) holds, but that between any two consecutive 



collisions, the carriers undergo to a free-evolution described by a (possibly time- 
depedent) Hamiltonian Hs(t) := J2m^Smi't) which are local (i.e. no direct interactions 
between the carriers is allowed). Under this circumstance it is possible to show that 



Eq. (31) still holds in the proper interaction picture representation at the price of 



allowing the generators of the resulting master equation to be explicitly time dependent 



as in Eq. (35 ). 



To see this we first notice that under the assumption that the collision time At 



The need of removing the free evolution of S from the description of the system dynamics is 



clearly evident also in our case. Indeed the condition (30) is clearly incompatible with the presence 



of free local contributions in the Hamiltonians Hg e as they will correspond to terms of the form 



H 



{free) 



5 Ie, i.e. contributions A 
\r^))E = (X"(r?))B = 1. 



it) 
s 



with B^ being the identity operator which will yield 
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is much shorter than the time interval that elapses between two consecutive collisional 
events (i.e. At ^ r„ — x„_i), the unitary operator which describes the evolution of the 
m-th carrier in its interaction with £ is now given by 

f^iS := f^5,„,i?„Ks„(rn, r„„i) Us^,e„_, (36) 

where Us^^e^ are the collisional transformations, r„ is the time at which the n-th collision 
takes place, and where Vs^{Tn, t:„_i) := Texp[— i J^" ^ dt'hs^it')] is the unitary operator 
which describes the free-evolution of Sm between the (n — l)-th and the n-th collision 
(in this expression Texp[- ■ ■] indicates the time-ordered exponential which we insert to 
explicitly account for possibility that the hs^ will be time- dependent). Define hence the 
operators 

4:'^-^i(^-0)4iv^5„(r„,0), (37) 
and the Hamiltonian 

Hsn.,En ■■= VliT^,0) Hs^,EnVsArn,0) = Yl ^sf ® 4! , (38) 

e 

which describes the coupling between Sm and S in the interaction representation 
associated with the free evolution of Sm. Notice that the operators ' are explicit 



functions of the index n which labels the collisions as in the case of Eq. (33) (here 



however the terms operating on 8 are kept uniform). Observing that for all i one has 



Vs^{Te,n-i)Vs^{Ti_i,Te_2) = '^5„(t-£, r^-2) we can now write Eq. (36) as 



f^£^.:=V5jr„,0)f/£„ (39) 

where U^sIIe unitary that defines the collisions of Sm with the sub-environments 

in the interaction representation, i.e. 

^S^,E„ Usrn,E^-i ■ ■ ■ Usm,Ei, (40) 

with 

Us„.,E„ = exp[-ig Hs^,En^t\ . (41) 
Similarly we can express the super-operators VV*^"'*^^ as 

y^inM) = V5(r„,0) o , (42) 

■'S,E„ " ■ ■ ■ " 's,Ei ' 



W(--):=Cg£o...oC^g, (43) 



■■= Me, o Us,„e, o...oMe,o Us,,e, , (44) 

with V5(r„, 0) being the super-operator associated with the joint free unitary evolution 
obtained by combining all the local terms of the carriers, i.e. V5(r„,0) : = 
(r„, 0) ■ • ■ Vsj^j {jn, 0). Defining hence R{n) the state of S and of the first elements of 
£ after n collisions in the interaction representation induced by Ks(r„,0) as 

R{n) = Vj(r„, 0) R{n) Ks(rn, 0) , (45) 



Master equation for cascade quantum channels: a collisional approach 14 

we get a recursive expression analogous to Eq. (9) with C^e„ i replaced by C^^^^^, i.e. 

Rin + l)=ci';ijRin)^v). (46) 
More precisely this expression coincides with that which, as in the case described at the 



end of Sec. |2.1[ one would have obtained starting from a collisional model in which no 
free evolution of the carriers is allowed but the collisional events are not uniform. Indeed 



the generators of the dynamics Hs^,e„ do have the same form of the Hamiltonians (33) 



Following the same prescription given there, we can then get an expression for the 
reduced density matrix p{n) = {R{n))£ which represents the state of the carriers after n 
collisions in the interaction picture with respect to the free evolution generated by Hs{t). 



operators Ag^{t) : 



Enforcing the limit (29) under the condition (34), one can verify that p{t) obeys to a 

being replaced by the time-dependent 



ME analogous to Eq. (|31|) with the operators Ag 

in/) 



lim 



At-s>0+ ^5, 



3.2. Enforcing the stability condition via a global unitary mapping 

Now that we have learned how to deal with free local evolutions terms operating between 
the collisional events, we show how to use this result to effectively enforce the stability 



condition of Eq. (30 ) for models in which it doesn't apply rigorously. Specifically we shall 



see that such condition can be imposed by first moving into an interaction representation 
with respect to a rescaled local Hamiltonian for the system S which maps the problem 



into one equivalent to that discussed in Sec. |3]T 
Indeed let 



H: 



Sm 



B 



(47) 



be the Hamiltonian which describe the collisions between the carriers and the sub- 
environments (notice that we are allowing the operators A^g^^ to depend explicitly from 
the n label to account for possible local free evolution of the carriers as discussed in the 



previous section). Suppose then that Eq. (30) does not hold. In this case we define 



Be 



(m,£) 



(^) 



Be 
{Be 



iv))E 



and write, 



where 



H 



An) 



AH 



Q TP 



E 



m ^5, 



in/) 



is a local Hamiltonian on Sm while 



AH 



Q TP 



A 



{n,e) 



B 



{m,l) 

En 



(48) 
(49) 

(50) 

(51) 

(52) 



is a rescaled coupling Hamiltonian. Differently from the original one given in Eq. (47), 
but similarly to Eq. ^33), it is built from operators -B^'^'' which explicitly depend 
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on the label m of the carrier 5*^, and which by construction satisfy the generalized 
condition (34), i.e. 



(53) 

Passing then in the interaction representation with respect to we can thus express 



the unitary evolution induced by Hs^^e^ 

, J /-At 

Us^,E,. = eM-wHsr^,E^^t] = e-'^^'sr.^' T eM'W / dt' AHs^,EAt')] , 

^0 

— ■ h A/ 

where e~^^ is a local unitary on Sm while 
Therefore the rhs of Eq. ([T]) can be written now as 

(54) 



where and Us^^Ej are the following unitary operators 

l(2) 



T/(«) . At -igh'^"''^'' At ^-igh''? At ^-igh'^}'' At 



and 



At 



Texp[-^g I dt' AHs^,EAt')]] ^t'^ ■ 



(55) 



(56) 



For future reference it is worth anticipating that the term within the round brackets 
admits the following expansion in At, 

Is^,E^ - z{gAt) AHs^,E, - ligAty {AHs^,e,)' - \{gAtf Qs„^,e, + 0{At') , 

where the last contribution originates from the time-ordering in the exponential of 
Eq. (56) and it is defined in terms of the first derivative of AH s^^E^it), i.e. 



dAHs^,EAt) 



dt 



t=o 



h'^^\AHs E 



(57) 



This yields the following expansion for the super-operator Us^,Ej associated to the 
unitary tJs^.Ep 

Us^,E, = Xs^,E, + {g^t) U's^^E^ + {gAtf Q'^^^E, + 0{{gAtf) , (58) 
with Ts^^E, being the identity map and with 



— I 



AH. 



i('L.EX---) ■■=AHs.A---)AHs.^ 



Am 



QSm,E., 



(59) 
, (60) 



where for easy of notation in this expression given a generic operator Qs,n,Ej we used 
the notation Qsm,Ej to represent its evolution via the unitary vj^ ^\ i.e. 



(61) 



Master equation for cascade quantum channels: a collisional approach 



16 



The expression Eq. (60) should be compared with Eq. (15): we notice that due to the 



presence of the time-ordering in Eq. (56) an extra term is present in the decomposition. 



We shall see however that when tracing out the sub-environments, such term plays no 



role in the evolution of the carries (see Eq. (75) below). 

With the above identities the row super-operator Tl^sl £ entering in Eq. (5) can thus 
be expressed as 



(62) 



where as usual V^"^ and Us^,Ej represent the super-operators associated with the unitary 

in) ~ 

transformations VgJ^ and Us^,Ej respectively, and where 

ns^,e ■■= -M®" o Us^^E^ 0...0 Us^,E. . (63) 
Accordingly Eq. ^ becomes 

(64) 



yy(n,Af) _ y(n) ^ yy{n,M) 
o 



with Vc"'' being the super-operator associated with the joint unitary 



(n) 



and with 



Sm 

:= iis^^ ^ o ns,,_,,s o ■ ■ ■ o iZs^^s o TIs^e . (65) 

This can also be written in terms of column super-operators as in Eq. ([T]). In particular 
we get 

= Cs,E^ o Cs,E^_, o ■ ■ ■ o Cs^E, o Cs,E, , (66) 

with 

Cs,E, ■■= Me, o I^Sm,Ej o ■ ■ ■ o Me, o I^SuE, ■ (67) 

Defining now 

R{n) := VP{R{n)) = [V^""^]^ R{n)vj''^ , (68) 

the state of the carriers and of the first n sub-environemental state in the interaction 
picture representation induced by the unitary V^^"'^ we have 

R(n + 1) = W("+i'^^)(p(0) ® r/®"+i) = Cs,E„^,{Rin) ® 77) . 



(69) 



Take then the partial trace over S of this expression and use Eq. (58) to expand Cs,Er,+i- 
Defining p{n) = {R{n))g we get 

(70) 



p{n + 1) = -p{n) + ((7 At) (C^,^„^, (i?(n) ® r]) 



+ {g/\tf (C^,^„,,(i?(n) ® v))^ + OiigAty 



where C[ 



S,E„ 



and a 



S,E, 



n+l 



are respectively the first and second order term of the 



expansion of Cs,e„+i, i-e. 



M 



C' 



m=l 



M-m+l 



E, °^Sm,E,°ME 



m—1 



^S,E, ■— 'S,Ei "T '-'S,Ei ' 



(71) 
(72) 
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M 



M-m+1 ^ 7 /// 



M M-l 



c 



n,b 
S,Ei 



oU'a p O 



m—1 



E E M%-'^' o U's^^^^^ o M%- o U',^^^^ o M^-' . (73) 



m'=m+l m=l 



As in the case of (11) one can verify that the first order contribution nulhfies. Indeed 
we have 



= 0, 



(74) 



m i 



because of Eq. (53). The remaining terms can be computed as in Eq. (21) and (24). 



Here we only stress on the fact that the component of Wg^ ^. that depends upon the 



operator Qsm,Ej (i-e. the extra term of Eq. (60)) do not contribute in the final result 



Indeed they only enters in the definition of Cc'*^ and produce the following term 



2 



(75) 



f Ed l^l^^^^^S^.iJn + l 



m £ 



m £ 



(n+1) Un+1,£) 



i?(n)®A^--\(r/)^ 

,i?(n)®A^--\(r/) 
rp{n) (Ol-(r/) 



. 



In summary Eq. (jTOj) yields 

p{n + 1) - p{n) g^At 



At 



1 



Y^C^{-p{n))+ E ^i"i'(P») + 0(^^At^) , (76) 



m'>m 



where now 
£ (■ 



£1' 



2A[ 



_^(n+l/)^(n+l/')/ 



and 



in,m' > 



{£,£') 



/ J im,m' 



£,£' 



x^(n+l/)^(n+l/) 

/ Sm Sin 



(■■■), At'-" 



A 
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In these expressions the coefficients 7™' and Tm^, differ from those in Eqs. (23) and 



(26) and are expressed by 

7^:2 := 7 ( Bf^'^M-'- ( Bt^'^ M-\ri) ) ) ^ . (77) 



E 



Also the operators A^"'^^'^ ■* stands for the operator A^^^'^ ■* in the interaction 



representation (68) induced by the transformation V^^"'^ i.e. 

:= [Vj")]t . (78) 

This follows from the fact that according to our definitions 

^s^.E^,. = [vt^]^&s^,E„,,vtj = [V^^]^ Qs^,E^,, V^i"^ =: Qs^,E^,. ■ (79) 



Taking now the limit (29) this finally yields a differential equation for p{t) as in 
Eq. (31) with operators being replaced by A^^(t) := limAt-s>o+ ^5^^'^^- It is worth 
noticing that in the continuos limit the transformation V^^"^ which define the interaction 
representation becomes: 



= ^mTexp[-ig / hs^{t')dt'] 
Jo 

with hs^{t) := limAi-^o+ h^sl- 



(80) 



4. Conclusions 

In the present manuscript we have reviewed some of technical aspects of the new 
method recently introduced in Ref. [T7] which allows one to derive in a consistent way, 
general master equation for cascade quantum system (i.e. multipartite quantum system 
which are unidirectionally coupled via a partially incoherent mediator). In particular 



we focused on the main assumption of the model (the stability condition of Eq. (30)) 
showing that it can be lifted by properly moving into a interaction picture representation 
with respect to the free dynamics of the system of interest. 
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